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ABSTRACT 



Comparison of the sidereal angular velocity of subphotospheric 
layers and small bright coronal structures during the declining 
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^\ ' Context. We compare solar differential rotation of subphotospheric layers derived from local helioseismology analysis of GONG++ 

dopplergrams and the one derived from tracing small bright coronal structures (SBCS) using EIT/SOHO images for the period August 

O !■ 2001 - December 2006, which correspond to the declining phase of solar cycle 23. 

Aims. The study aims to find a relationship between the rotation of the SBCS and the subphotospheric angular velocity. The north- 
south asymmetries of both rotation velocity measurements are also investigated. 

Methods. Subphotospheric differential rotation was derived using ring-diagram analysis of GONG++ full-disk dopplergrams of 1 min 
cadence. The coronal rotation was derived by using an automatic method to identify and track the small bright coronal structures in 
^ i ' EIT full-disk images of 6 hours cadence. 

Results. We find that the SBCS rotate faster than the considered upper subphotospheric layer (3Mm) by about 0.5 deg/day at the 
equator. This result joins the results of several other magnetic features (sunspots, plages, faculae, etc.) with a higher rotation than the 
solar plasma. The rotation rate latitudinal gradients of the SBCS and the subphotospheric layers are very similar. The SBCS motion 

■ shows an acceleration of about 0.005°day "'/month during the declining phase of solar cycle 23, whereas the angular velocity of 
' subsurface layers does not display any evident variation with time, except for the well known torsional oscillation pattern. Finally, both 

subphotospheric and coronal rotations of the southern hemisphere are predominantly larger than those of the northern hemisphere, 
^-j- . At latitudes where the north-south asymmetry of the angular velocity increases (decreases) with activity for the SBCS, it decreases 

• ■ (increases) for subphotospheric layers. 

■ 
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1 . Introduction increases with depth in the upper layers of the convection zone, 

at intermediate latitudes. 

The solar angular rotation velocity is a function of latitude, time, The coronal angular velocity has been measured using sev- 
53 1 and height above or depth below the solar photosphere. This phe- eral coronal structures such as coronal green lines (e.g, Altrock, 
' nomenon is known as the solar differential rotation. The most 2003; Badalyan et al. 2006), radio emission flux (e.g, Mouradian 
commonly established latitudinal dependence of the angular ve- et al., 2002), and coronal holes (e.g, Insley et al., 1995). Many 
locity is given by the empirical relation 0.(0) - A + Bsin 2 + of these features show two rotational modes where high-latitude 
Csin 4 6», where 9 is the latitude, A is the equatorial rotation, B regions rotate more rigidly than low-latitude regions. However, 
and C are the differential rotation coefficients, and C is usually P art of this behaviour has been related to the lifetime of the 
neglected for low and mid-latitude measurements. In terms of structure and not specifically to that of the corona (see a re- 
amplitude, the rotation of the quiet Sun regions, given for in- view b Y Schroter, 1985). Moreover, small bright coronal struc- 
stance by Doppler-shift measurements of Fraunhofer lines, is tures (SBCS, hereafter) have been used to estimate the coro- 
found to be slower than the rotation of the photospheric mag- nal rotation from data recorded by the EUV imaging telescope 
netic tracers (individual sunspots, sunspot groups, faculae, su- (EIT) on board the s °l ar and heliospheric observatory (SOHO). 
pergranules, etc.). Moreover, magnetic features exhibit different These features are mostly short-lived structures with a lifetime 
angular velocities depending on their evolution, age, and size up to 54 hours (Brajsa et al., 2008) and are formed at a height 
(see a review by Beck, 2000). So far, this difference is explained of about 8000-12000 km above the photosphere (Brajsa et al., 
by magnetic features being rooted at different depths (Ruzdjak 2004). Brajsa et al. (2001,2002) explored these structures using 
et al., 2004). Several local and global helioseismology analyses 284 A EIT filtergrams with a 6 hr cadence, whereas Karachik et 
of solar acoustic modes, such as time-distance helioseismology al. (2006) used 194 A EIT filtergrams with 12 min cadence. Both 
analysis (e.g., D'Silva, 1996) and f-mode analysis (e.g., Corbard authors confirm the differential rotation of the corona through 
& Thompson, 2002), have confirmed that the angular velocity the SBCS with a close relation to the differential rotation of at- 
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Table 1. Unsigned differential rotation parameter B derived from the linear fit of the rotation averaged over 3 time frames [Aug 2001 
- Dec 2006; Aug 2001 - Mar 2004; Apr 2004 - Dec 2006] and sin 2 (6») for all latitudes [-45°,45°], southern hemisphere [-45°,0°] 
and northern hemisphere [0°,45°]. 



August 2001 -December 2006 



August 2001 -March 2004 



April 2004-December 2006 





All latitudes 


North 


South 


All latitudes 


North 


South 


All latitudes 


North 


South 


SBCS 


2.72 


±0.04 


2.77 ± 0.04 


2.64 


±0.04 


2.70 ± 0.03 


2.70 ± 0.04 


2.69 ± 0.04 


2.73 


±0.07 


2.82 


±0.07 


2.59 


±0.08 


3 Mm 


2.77 


±0.01 


2.74 ± 0.01 


2.73 


±0.01 


2.77 ± 0.01 


2.77 ± 0.02 


2.70 ± 0.02 


2.77 


±0.01 


2.72 


±0.02 


2.76 


±0.02 


6 Mm 


2.78 


±0.01 


2.77 ±0.01 


2.73 


±0.01 


2.78 ± 0.01 


2.80 ± 0.02 


2.70 ± 0.02 


2.78 


±0.01 


2.74 


±0.02 


2.76 


±0.02 


15 Mm 


2.88 


±0.01 


2.90 ± 0.01 


2.80 


±0.01 


2.88 ±0.01 


2.93 ± 0.02 


2.77 ± 0.02 


2.88 


±0.01 


2.87 


±0.02 


2.82 


±0.02 



mospheric magnetic features. Moreover, the SBCS are interest- 
ing features to be used for the rotation estimation along the cycle 
since they appear at both minimum and maximum cycle phases, 
in contrast to other activity features that are absent at minimum 
activity or, for instance, coronal holes with their polar concen- 
tration at that period. Their relation to the photospheric mag- 
netic field has been studied by Pres and Phillips (1999), who 
find that the time evolution of bright coronal points observed 
with EIT/SOHO is correlated with the rise and fall of the mag- 
netic field given by the MDI/SOHO magnetograms. 

In this paper, we compare the rotation velocity of subpho- 
tospheric layers (from 3Mm to 15Mm) obtained from ring- 
diagram analysis of GONG++ data using the GONG ring- 
diagram analysis pipelin43 (Corbard et al., 2003) with the one 
measured by tracing the SBCS observed in EIT/SOHO images 
(at 284 A) using the automatic method described in Brajsa et al. 
(2001). The purpose of this comparison is to get a complemen- 
tary view of the relation between coronal features and their root 
layers below the photosphere. Moreover, we investigate the lati- 
tudinal gradient of the angular velocity to check the assumption 
of the two rotational modes and the variation of the latitudinal 
profile of the coronal rotation during the declining phase of the 
solar cycle 23 as reported by Altrock (2003). Finally, a com- 
parison between the rotations of the northern and the southern 
hemispheres is shown for both SBCS and subphotospheric lay- 
ers. 



2. Data reduction 

2.1. Rotation of subphotospheric layers 

Subphotospheric angular velocity is measured using ring- 
diagram analysis which is a local helioseismology technique 
based on frequency-shift measurements of high-degree acoustic- 
modes to infer horizontal velocity flows at different subphoto- 
spheric depths (Hill, 1988). By dividing the full-disk image into 
189 overlapping patches with centres separated by 7.5° in lat- 
itude and longitude with a latitude range of [-52.5,52.5] and a 
central meridian distance range of the same extent, each region 
is remapped and tracked at the surface rotation rate of its central 
latitude in consecutive images with a 1 min cadence. The result- 
ing data cubes are Fourier-transformed, and the resulting power 
spectra are fitted with a Lorentzian profile to derive frequency 
shifts that are directly related to horizontal velocity components 
by assuming a plane wave approximation. The horizontal veloc- 
ity as function of depth is then deduced using inversion methods. 



http://gong2.nso. edu/archive/patch.pl?menutype=h 



We considered the depth range [3Mm,15Mm]. We avoided 
extending the radial depth to upper layers since surface measure- 
ments of the flow are affected by poor estimation of the inver- 
sion kernels because of the complexity of the photospheric lay- 
ers. Deeper measurements were not possible because the small 
size of the patches only gives access to high spatial frequencies 
(high spherical harmonic degree). The sidereal angular velocity 
for each local area (latitude, Carrington longitude) was recov- 
ered by adding the measured east-west flow (zonal flow) to the 
tracking rate. Finally, the rotation rate at each latitude is given by 
averaging over all the available Carrington longitudes for each 
month. The synodic-sidereal correction was taken into account 
for each day. 

2.2. Rotation of small bright coronal structures 

The data reduction to determine the motions of bright coro- 
nal structures was described in detail by Brajsa et al. (2001). 
Although not just coronal bright points (CBP) were selected in 
the cited paper, this abbreviation was used. In fact, besides points 
also small loops and small active regions were selected. We pre- 
fer to call the used tracers small bright coronal structures (SBCS) 
now. The interactive method has up to now only been applied to 
images obtained in 1998 and 1999, because it is very time con- 
suming. An automatic method was developed later to cover more 
images from most of the years of an activity cycle. The main idea 
was to automatically select small structures in three consecutive 
images and find whether they belong to a stable structure within 
12 hours so they only show small changes in latitude and in lon- 
gitude. Slight differences were applied to the parameters of the 
automatic reduction program (Brajsa et al. 2001). First, the cir- 
cumferences of the allowed structures was changed to the range 
of 30 to 80 pixels, which equals 80 to 210 arcseconds. The range 
of relative intensities was 100 to 600 units. Most important was 
a change in the allowed difference in latitude from one image 
to the next from 1 degree to 0.8 degree and the application of 
a differential rotation for the allowed motions in longitude: The 
limits were 9.5 to 16.5 degrees per day of the synodic rotation, 
but these limits were reduced by -3.0*sin 2 (#), where 9 is the 
latitude. 

The conversion of synodic rotation rates to sidereal was per- 
formed for each observing time separately. Filtering after a first 
fit to the differential rotation law with C — was applied in the 
same way as described in Section 4.2 of Brajsa et al. (2001). This 
filtering was applied only when data from certain periods (e.g. 
months, years or longer periods) are combined. Thus the sum 
of all structures used is slightly different when adding monthly 
numbers and comparing with annual sums. At the beginning the 
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Fig. 1. Rotation velocity derived from ring-diagram analysis at 
3Mm (asterisk), 7Mm (cross), 15Mm (diamond) (lines from 
dark to light grey) and SBCS tracing (black) averaged over the 
periods: August 2001 - March 2004 (dashed), April 2004 - 
December 2006 (dashed dotted). 

data were selected for all latitudes in each month and only differ- 
ential rotation parameters fitted using C — 0. To adapt the reduc- 
tion to the GONG++ data, a new program was added, where the 
results were sampled for overlapping bins of 15 degrees width 
in latitude and a stepping of 7.5 degrees. Thus a data cube was 
derived with rotation rate results and their errors for 15 latitude 
bins and 65 months from August 2001 until December 2006. 

3. Solar differential rotation 

We measured the subphotospheric and coronal rotation veloci- 
ties for the period August 2001 - December 2006, which corre- 
sponds to the declining phase of solar cycle 23. GONG++ data 
covers all the months with a good duty cycle. The variation in the 
zonal flow (i.e rotation residual) as derived from ring-diagram 
analysis of GONG++ data during this period has been studied 
in detail by Komm et al. (2009) with 3 extra months in 2007. 
They find a positive correlation between the unsigned magnetic 
flux and the zonal flow and report a higher zonal flow in patches 
with activity than for the quiet Sun patches. They do not find any 
particular long-term pattern of the zonal flow except the well- 
known pattern of torsional oscillations (Howard and LaBonte, 
1980). However, a slight increase in the zonal flow in patches 
with activity was observed at minimum activity. Figure Q] shows 
the sidereal angular velocity given by the coronal and subphto- 
spheric measurements for the latitudinal range [-52.5°, 52. 5°]. 
From ring-diagram analysis measurements, we see that deep 
subsurface layers rotate faster than those closer to the surface 
with a decreasing radial gradient as the latitude increases. The 
radial gradient variation in the rotation velocity estimated from 
ring-diagram analysis of the same data was studied in detail by 
Zaatri and Corbard (2009), who report that the use of higher lati- 
tudes can lead to a reversal of the radial gradient from negative to 
positive. Moreover, Fig. Q] shows a smaller northern radial gra- 
dient than the southern gradient at high latitudes. From SBCS 
tracing, the rotation rate is found to be higher than for all the 
considered subsurface layers by about 0.5°day~' at the equator 
(this corresponds to a velocity of about 70 m/s, and a frequency 
of about 16 nHz) compared to that of the subsurface layer at 
3Mm depth. This difference is in good agreement with the angu- 




-0.005 



-40 -20 20 

latitude [deg] 



40 



Fig. 2. Upper panel: Slope of the linear fit of the rotation with 
time derived from ring analysis at 3Mm, 6Mm and 15Mm (from 
dark to light grey) and SBCS tracing (black) for the period 
Aug 2001 - Dec 2006. Lower panel: Slope of the linear fit 
of the SBCS angular velocity with time for the periods: Aug 
2001 - Mar 2004 (dashed line, triangles); Apr 2004 - Dec 2006 
(dashed-dotted line, squares). Squares are shifted by 1° to avoid 
overlapping of error bars. 



lar velocity difference between the photospheric active regions 
and their surrounding quiet solar plasma (e.g, Koch, 1984). 

The upper panel of Fig. [2] shows the angular acceleration 
of the SBCS and of the subsurface layers. The subphotospheric 
layers have no significant acceleration compared to that of the 
SBCS. However, it is worth mentioning that the absence of ac- 
celeration at low latitudes and the week acceleration at latitudes 
higher than 20° (up to 10~ 3 deg/day/month) is consistent with the 
torsional oscillation pattern for the years 2001 - 2006 (see Fig. 
30 of Howe 2009). The SBCS rotation shows a positive increase 
with an acceleration of about 0.005°day~ 1 /month. This accelera- 
tion also depends on activity as shown by the lower panel of Fig. 
|2] At higher activity, the acceleration is clearly seen to be lower 
than the lower activity for the southern hemisphere; however, it 
shows a more complicated latitudinal and activity variation in 
the northern hemisphere. It is worth mentioning that the photo- 
spheric magnetic features are also rotating faster at low activ- 
ity, as seen from several observations (e.g, Brajsa et al., 2006). 
Lastly, we investigated the latitudinal gradient of the rotation 
by measuring the differential rotation coefficient B. We avoided 
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taking the highest latitude ±52.5° where the SBCS are almost 
absent at low activity, leading to angular velocities which are 
mostly interpolated values. Table 1 gives the B values for both 
SBCS and three subphotospheric layers, and B is seen to increase 
with depth at the outer part of the convection zone. Moreover, 
neither the small SBCS nor the subphotospheric layers show any 
significant variation with the activity of the latitudinal gradient 
of their angular velocity. 

4. North-south asymmetry 

We investigated the north-south asymmetry of the coronal ro- 
tation velocity through the SBCS angular velocity change be- 
tween the northern and southern hemispheres by evaluating the 
difference between two symmetric bins. For instance, a rotation 
velocity averaged over a bin centred at 15° is subtracted from 
that derived from a bin centred at -15°, the difference value is 
shown at latitude 15°. This quantity is evaluated for the unsigned 
latitude range [7.5°,45°] in the upper panel of Fig. [3] The north- 
south asymmetry of the coronal rotation reveals a faster south- 
ern hemisphere during the period under analysis. Also shown 
is a latitudinal dependence of this property with an increasing 
north-south difference in rotation from the equator to about 20°, 
where the amplitude is more pronounced with strong activity. At 
latitudes higher than 20°, this asymmetry continues to increase 
at the low activity level (end of cycle 23), whereas it is more 
irregular at the beginning of the declining phase and becomes 
even positive around 30°. To see whether the asymmetry in the 
coronal rotation velocity is related to the coronal activity, we in- 
vestigated the activity asymmetry via the asymmetry index given 
by (N - S)/(N + S) for each pair of symmetric bins where N is 
the number of structures in the northern bin and S is the num- 
ber of structures in the southern bin. The asymmetry index and 
its statistical mean error are evaluated on a monthly basis and 
shown in Fig. [3] The figure shows that the northern coronal ac- 
tivity prevails during the first half of the declining phase of solar 
cycle 23 and the southern activity is more dominant during the 
second half of this period. This agrees with the flare activity dur- 
ing this cycle, which favoured the north at maximum activity and 
the south during the declining phase (Joshi et al. 2007). Figure[3] 
clearly shows that, in the first half of the declining phase of solar 
cycle 23, the northern hemisphere, which is more active, rotates 
more slowly than the southern hemisphere, and it remains slower 
at the end of the cycle even though it becomes less active than 
the southern hemisphere. This disagrees with several studies re- 
porting that the active hemisphere rotates more slowly than the 
less active one (Obridko et al. 2001). 

Figure [4] shows the difference between the northern and 
southern rotation velocities for three subsurface layers at depths 
of 3Mm, 7Mm, and 15Mm. As for the SBCS, subsurface rota- 
tion velocity is predominantly southern except at the intermedi- 
ate latitudes of the upper subsurface layers during the low activ- 
ity epoch. This southern dominance also characterised the pho- 
tospheric magnetic activity during solar cycle 23 as shown from 
sunspot indicators (See Li et al. 2008 and references therein). 
Besides, Zaatri et al. (2006) have analysed the first three years of 
our GONG++ dataset and found that the zonal flow pattern and 
its north-south asymmetry are strongly related to the distribution 
of magnetic activity. Moreover, a clear increase of this asym- 
metry with depth is seen for the whole period. This confirms, 
for a longer period, the results by Zaatri et al. (2006) who found 
the same north-south asymmetry depth dependence for the zonal 
flow. Unlike the SBCS, the north-south asymmetry in the sub- 
surface rotation velocity shows a slight increase with decreasing 
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Fig. 3. Upper panel: Difference between the northern and south- 
ern angular velocity of SBCS averaged over the periods August 
2001 - December 2006 (full lines), August 2001 - March 
2004 (dashed,triangles), April 2004 - December 2006 (dashed- 
dotted,squares). Lower panel: Asymmetry index of the coronal 
activity derived from the number of observed SBCS in the north 
(N) and in the south (S) given by the formula (N-S)/(N+S) by 
summing the structures over the periods given in the upper panel. 
The thin grey line shows the zero level. Triangles and squares are 
shifted by 1 to the left and right to avoid overlapping the error 
bars. 

activity, at low latitudes. At latitudes higher than 20° and close 
to solar cycle maximum, the subsurface asymmetry is slightly 
changing until around 40°. During the low activity epoch, the 
rotation velocity of the southern hemisphere gets very close to 
the northern velocity and even smaller around latitude 30° for 
the upper subsurface layers. 

Some hints of the north-south asymmetry are also given from 
the latitudinal gradient in Table 1 where the subsurface layers 
show more varying B values with depth in the northern hemi- 
sphere than in the southern hemisphere for both activity levels. 
The SBCS have a more pronounced north-south variation of B 
at low activity than at high activity. 

5. Discussion and conclusion 

In this work, we compared the rotation of the corona as revealed 
by tracing its most profuse features, the SBCS, and that of the 
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angular velocity of subsurface layers as a function of unsigned 
latitude and for depths 3Mm, 7Mm and 15Mm. Shown are the 
periods : August 2001 - December 2006 (full lines), August 
2001 - March 2004 (dashed), April 2004 - December 2006 
(dashed-dotted). The thin grey line shows the zero level. 



uppermost layers of the convection zone. We report the follow- 
ing results 

(i) There is a difference of about 0.5°day _1 between the an- 
gular velocity of the SBCS and the upper considered subpho- 
tospheric layer at 3Mm below the surface. This difference de- 



creases while considering deeper layers due to the increase of 
the angular velocity with depth at intermediate latitudes. At the 
highest latitudes of our considered range [-52.5°, 52.5°], the ra- 
dial gradient of the rotation velocity decreases considerably and 
the difference between the angular velocities of the SBCS and 
the various subsurface layers becomes very similar. 

(ii) The temporal variation of the SBCS angular velocity dur- 
ing the declining phase of the solar cycle 23 shows an acceler- 
ation of about 0.005°day "'/month on average and reaches its 
maximum during the low activity period (end of cycle 23) in the 
southern hemisphere. It is not clear whether this acceleration is 
related to the chosen tracers or if it is the property of the corona 
and whether it is specific to the studied period. So far, there have 
been no studies of the coronal rotation using other activity fea- 
tures for the last years of solar cycle 23. However, results from 
other periods have reported a deceleration that makes the coro- 
nal rotation tend to that of the photosphere at solar minimum (e.g 
Giordano et al. 2008 for the period 1996-1997), whereas Mehta 
(2005) did not find any evidence of the coronal rotation varia- 
tion during four consecutive solar cycles (19 to 22). The angular 
velocity does not show a significant temporal variation except 
a slight acceleration at high latitudes which corresponds to the 
torsional oscillation signal during the period 2001-2006. 

(iii) The latitudinal dependence of the angular velocity from 
the two measurements was studied through the differential rota- 
tion parameter B as given by the most commonly used rotation 
law at intermediate latitudes (Q(0) = A + Bsin 2 6). We found 
a very similar latitudinal gradient B of coronal and subphoto- 
spheric rotations, except that the unsigned value of B slightly 
increases with depth. From this point of view, we do not see 
any particularity in the latitudinal behaviour of the coronal ro- 
tation as mentioned by several authors from long-term coronal 
features observations, notably, the widely mentioned existence 
of two different rotational modes with a more rigid rotation of the 
corona than for the photosphere at low latitudes (Altrock 2003). 
Clearly, the SBCS behave more like magnetic photospheric fea- 
tures which endorses the hypothesis of the formation of EUV 
bright points as being caused by the heating derived from the 
dissipation of electric currents that are formed in the solar at- 
mosphere from the displacement of footpoints of magnetic flux 
tubes by the photospheric plasma motion (see Santos et al. 2008 
and references therein). 

(iv) For both corona and subphotospheric layers, the south- 
ern hemisphere rotates predominantly faster than the northern 
hemisphere. The north-south velocity difference increases with 
depth. This result was already shown in Zaatri et al. (2006) and 
is now confirmed for a longer period including the very low ac- 
tivity period of cycle 23. Also, when the north-south asymmetry 
of the rotation velocity is lower than its temporal average for 
the subsurface layers, it is higher for the SBCS and vice-versa. 
Clearly, the southern coronal hemisphere rotates faster than the 
northern hemisphere regardless of which hemisphere is more ac- 
tive. It would be interesting to see whether this persists on a time 
scale of multiple cycles as opposed to the relatively short time 
scale of our study. 
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